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' Abstract 

o . 

, The low-laying frequencies of characteristic quasi-normal modes (QNM) of Schwarzschild-de 

Sitter (SdS) black holes have been calculated for fields of different spin using the sixth-order 
^ ' WKB approximation and the approximation by the Poshl- Teller potential. The well-known 

asymptotic formula for large I is generalized here on a case of the Schwarzschild-de Sitter black 
hole. In the limit of the near extreme A term the results given by both methods are in a very 
good agreement, and in this limit fields of different spin decay with the same rate. 
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1 Introduction 



Within general relativity it is well-known that there are three stages to the dynamical evolution of 
field perturbation on a black hole background: the initial outburst from the source of perturbation, 
the damping (quasi-normal) oscillations and asymptotic tails at very late time. The evolution signif- 
icantly depends on the asymptotic behavior of the background: for different values of cosmological 
constant A there are three types of behavior depending on the sign of A. Thus, the remarkable fact 
is that while for asymptotically flat backgrounds we have inverse power law tails, for asymptotically 
Anti-de Sitter background it is exponential low that governs the decay of a field at all times. The 
same exponential decay at very late times was observed by Brady et al [1] for asymptotically de 
Sitter background. 

Apart from original interest coming from possibility to observe quasi-normal ringing with the help 
of gravitational wave detectors (see, e.g., [2] for review), the quasi-normal modes bring now a lot of 
interest in different contexts: in Anti-de Sitter(de Sitter) /Conformal Field Theory correspondence 
[3, 4, 5], when considering thermodynamic properties of black holes in loop quantum gravity [6], in 
the context of possible connection with critical collapse [3, 7]. 

The gravitational QNMs of Schwarzschild-de Sitter black hole were considered in [8, 9, 5]. Re- 
cently, asymptotic values of QNMs were studied for nearly extremal Schwarzschild-de Sitter black 
hole and it was shown that in the near extremal limit the effective potential of the wave equation 
reduces to the Poshl- Teller potential, and it is expected therefore the Poshl-Teller method provides 
the best accuracy near the extremal value of A [11]. 

QNMs corresponding to decay of fields of different spin are extensively studied for Schwarzschild 
and Schwarzschild-anti de Sitter backgrounds. At the same time there is lack of such study of 
the Schwarzschild-de Sitter case. This motivated us to calculate QNMs of scalar, electromagnetic, 
gravitational and Dirac fields for Schwarzschild-de Sitter black hole. We also compare the results 
found here with those obtained recently in [10, 11]. 

In section 2 we give the perturbation equations to describe fields of different spin in SdS back- 
ground. Section 3 presents comparison of QNM frequencies obtained through sixth order WKB 
formula and approximation by the Poshl-Teller potential. 

There was shown that the frequencies of the massless Dirac field perturbations [10] are the same 
for both chiralities in SdS background as well. Near the extremal value of A the sixth order WKB 
formula gives the values which are very close to those obtained within the Poshl-Teller potential 
approximation. This agrees with recent work of Cardoso and Lemos [11]. 

2 Perturbation equations and methods 

It is well known that the Schwarzschild-de Sitter black hole of mass M is described by the metric: 



A perturbation equation on this background can be reduced to the Schrodinger wave-like equations: 



ds' = f{r)dt'' 



dr'^ 




fir) 




(2) 
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where we assume exponential dependence of ^ on time (~ e "^*), u — LORe — icoim, r* is the "tortoise" 
coordinate: 

- m (3) 

Under the positive real part QNMs, by definition, satisfy the following boundary conditions: 

*(r*) ~ C±exp(±ia;r*), r — > ±00, (4) 

which correspond to purely out-going waves at infinity and in-going waves at the horizon. 
The corresponding perturbation equations can be rewritten in form (2) with potentials 

V. - /M^; (6) 

2f{r)9M^ + 3c'^Mr^ + c^{l + cy + 3M^{3cr -Ar^) _/(/ + !) 

Veven - {3M + Cr)^ ' ~ ^ ' ^ ^ 

V, . ^J^M^r) ± - D); (9) 

for massless scalar, electromagnetic, gravitational (odd end even parities) and Dirac fields respectively 
(see, e.g., [12] and references therein), where I is the multipole number, k — I + 1 ior (+) sign and 
K — I ior (-) sign in (9). All the potentials vanish at the horizon and tend to a constant at infinity, 
having a peak near the horizon. It allows to use the Poshl-Teller potential approximation as well as 
WKB method. 

The WKB approach is based on the analogy with the problem of waves scattering near the peak 
of the potential barrier V{r*) in quantum mechanics, where uj^ plays a role of energy. The approach 
was used by Schutz and Will [13], developed to the third order by Iyer and Will [14], and recently 
developed by Konoplya to the sixth order beyond the eikonal approximation [15]. The result has the 
form: ^ 

i^^^^^^^-A2-A3-A4-A5-A6 = n+ J, (10) 



where Vq is the height and Vq is the second derivative with respect to the tortoise coordinate of the 
potential at the maximum. A2 and A3 can be found in [14], A4, A5 and Ag are presented in [15]; the 
corrections depend on the value of the potential and higher derivatives of it at the maximum. 

Generally, WKB method was effectively used for finding quasi-normal frequencies in a lot of papers 
[16]. Comparing the values calculated by the difi^erent order WKB formula one can judge about the 
convergence to the some unknown accurate result as the WKB order increasing. We expect that the 
possible error is less then the difference between the sixth and the fifth order WKB values^. Note, 
that standard machine precision in MATHEMATICA gives large errors due to numerical truncations. 
To exclude sueli " noises" i^ e liad to hold 40 digits in all intermediate values. 

^The accuracy of values are to take as a half of the number last order (e. g. 0.125 means 0.125 ± 0.003), but it is 
much better for most of them. 
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For comparison we will use the semi-analytical method for calculations proposed by B. Mashhoon 
[18] who used Poshl- Teller approximate potential 



PT 



cosh^(r*/6) 



(11) 



It contains two free parameters (Vq and b) which are used to fit the height and the second derivative 
of the potential V{r*) at the maximum [18]. 

The quasi-normal modes of the Poshl-Teller potential can be evaluated analytically: 




1 
4 



{n+-)t 



(12) 



3 QNMs of SdS black holes 



Both methods give best result for low overtone number but it is known that the higher modes give 
larger values of the imaginary part of frequencies [2]. Therefore these modes decay faster then low- 
lying ones that are more interesting for us. So we calculated some first quasi-normal mode frequencies 
for massless scalar, electromagnetic and gravitational and Dirac fields using both of the methods. 
The results are presented below (all values are measured in the black hole mass (M) units). 



QNMs of massless scalar field: I — 1, n — 



A 


^' (WKB) 


^' (P-T) 





0.2929 


- 0.0978^ 


0.299 


- O.IOH 


0.02 


0.2603 


-0.0911i 


0.263 


- 0.093i 


0.04 


0.2247 


- 0.0821i 


0.226 


- 0.083i 


0.06 


0.1854 


- O.OTOli 


0.187 


- omii 


0.08 


0.1404 


- 0.0542i 


0.141 


- 0.055i 


0.09 


0.11392 - 


- 0.04397i 


0.1147- 


- 0.0443i 


0.10 


0.08156 - 


' 0.03121i 


0.0819- 


- 0.0315i 


0.11 


0.02549 - 


- 0.00965i 


0.02550 - 


0.00967i 



QNMs of massless scalar field: I — 2, n — 



A 


^' (WKB) 


^' (P-T) 





0.48364 - 


0.09677^ 


0.487 


- 0.098i 


0.02 


0.43461 - 


0.08858i 


0.437 


- 0.089i 


0.04 


0.38078 - 


0.07876i 


0.382 


- 0.079i 


0.06 


0.32002 - 


0.06685Z 


0.321 


- 0.067i 


0.08 


0.24747 - 


0.05197i 


0.248 


- 0.052i 


0.09 


0.20296 - 


0.04256i 


0.2033 - 


- 0.0427i 


0.10 


0.14661 - 


0.03069i 


0.1468 - 


- 0.0308i 


0.10 


0.04617- 


0.00963i 


0.04617 - 


0.00963i 



QNMs of massless scalar field: I — 2, n — 1 



QNMs of electromagnetic field: I — 1, n — 



A 


u (WKB) 


u (P-T) 





0.46385 


- 0.29563i 


0.487 


- 0.294i 


0.02 


0.42084 


- 0.26862i 


0.437 


- 0.268i 


0.04 


0.37165 


- 0.23796i 


0.382 


- 0.238i 


0.06 


0.31422 


-0.20163i 


0.321 


- 0.202i 


0.08 


0.24426 


- 0.15627i 


0.248 


-0.156i 


0.09 


0.20098 


- 0.12794i 


0.2033 - 


- 0.1282i 


0.10 


0.14576 


- 0.09212i 


0.1468 - 


- 0.0923i 


0.11 


0.04614 


- 0.02889i 


0.04617 - 


0.02890i 



A 


u (WKB) 


u (P-T) 





0.2482 


- 0.0926i 


0.255 


- 0.096i 


0.02 


0.2259 


- 0.0842i 


0.231 


- 0.087i 


0.04 


0.2006 


- 0.0748i 


0.204 


- 0.077i 


0.06 


0.1709 


- 0.0639i 


0.173 


- 0.065i 


0.08 


0.1339 


- 0.0502i 


0.135 


-0.05H 


0.09 


0.11053- 


- 0.04156i 


0.1110- 


- 0.0419i 


0.10 


0.08035 - 


- 0.03028i 


0.0805 - 


- 0.0304i 


0.11 


0.02545 - 


- 0.00962i 


0.02546 - 


0.00962i 
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QNMs of electromagnetic field: I — 2, n — 



QNMs of electromagnetic field: I — 2, n — 1 



A 


uj (WKB) 


a; (P-T) 




A 


UJ (WKB) 


u (P-T) 







— u.uyouit 




- 0.096i 







0.43653 


— u.zyu / oi 




— u.zoyt 




0.41502 


- 0.08615i 


0.418 


- 0.087i 




0.02 


0.39900 


- 0.26202i 


0.418 


- 0.261i 


0.04 


0.36723 


- 0.07624i 


0.369 


- 0.077i 




0.04 


0.35602 


- 0.23065Z 


0.369 


- 0.231i 


0.06 


0.31182 


- 0.06478i 


0.313 


- 0.065i 




0.06 


0.30498 


-0.19516i 


0.313 


-0.196i 


0.08 


0.24365 


- 0.05067i 


0.244 


- 0.051i 




0.08 


0.24046 


- 0.15223i 


0.244 


- 0.153i 


0.09 


0.20085 


- 0.04180i 


0.2012- 


- 0.0419i 




0.09 


0.19907 


- 0.12549i 


0.2012 - 


- 0.1258i 


0.10 


0.14582 


- 0.03037i 


0.1459- 


- 0.0304i 




0.10 


0.14515 


- 0.09114i 


0.1459 - 


- 0.0913f 


0.11 


0.04614 


- 0.00962i 


0.04615 - 


0.00962i 




0.11 


0.04612 


- 0.02886i 


0.04615 - 


0.02887i 








QNMs of 


gravitational field: I — 2, n 


= 







A 


Lu (WKB) (odd) 


UJ (P-T) (odd) 


u (WKB) (even) 


UJ (P-T) (even) 




0.02 
0.04 
0.06 
0.08 
0.09 
0.10 
0.11 


0.3736 - 0.0889i 
0.3384 - 0.0817i 
0.2989 - 0.0733i 
0.2533- 0.0630i 
0.1975- 0.0499i 
0.16261 - 0.04136i 
0.11792 - 0.03021i 
0.03727 - 0.00962i 


0.378 - 0.091i 
0.342 - 0.083i 
0.301 - 0.074i 
0.254 - 0.064i 
0.198 - 0.050i 
0.1629 - 0.0415i 
0.1180 - 0.0303i 
0.03727 - 0.00962i 


0.3737 -0.0889i 
0.3384 - 0.0817i 
0.2989 - 0.0733i 
0.2533 - 0.0630i 
0.1975 - 0.0499i 
0.16261 - 0.04137i 
0.11792 - 0.03021i 
0.03727 - 0.00962i 


0.378 - 0.091i 
0.342 - 0.083i 
0.301 - 0.074i 
0.254 - 0.064i 
0.198 - 0.050i 
0.1629 - 0.0415i 
0.1180 - 0.0303i 
0.03727 - 0.00962i 


QNMs of gravitational field: / = 2, n = 1 


A 


UJ (WKB) (odd) 


UJ (P-T) (odd) 


UJ (WKB) (even) 


UJ (P-T) (even) 




0.02 
0.04 
0.06 
0.08 
0.09 
0.10 
0.11 


0.3463 - 0.2735i 
0.3185 - 0.2488i 
0.2858 - 0.2215i 
0.2457 - 0.1896i 
0.1941 - 0.1497i 
0.16079 - 0.12412i 
0.11724 - 0.09063i 
0.03725 - 0.02885i 


0.378 - 0.272i 
0.342 - 0.249i 
0.301 - 0.223i 
0.254 - 0.191i 
0.198 - 0.150i 
0.1629 - 0.1246i 
0.1180 - 0.0908i 
0.03727 - 0.02885i 


0.3467 - 0.2739i 
0.3188 - 0.2491i 
0.2858 - 0.2217i 
0.2457 - 0.1898i 
0.1941 - 0.1498i 
0.16079 - 0.12415^ 
0.11724- 0.09064i 
0.03725 - 0.02885i 


0.378 - 0.272i 
0.342 - 0.249i 
0.301 - 0.223i 
0.254 - 0.191i 
0.198 - 0.150i 
0.1629 - 0.1246i 
0.1180 - 0.0908i 
0.03727 - 0.02885i 



QNMs of gravitational field: I — 3, n — 



A 


UJ (WKB) (odd) 


UJ (P-T) (odd) 


UJ (WKB) (even) 


UJ (P-T) (even) 




0.02 
0.04 
0.06 
0.08 
0.09 
0.10 
0.11 


0.599443 - 0.092703i 

0.543115 - 0.084495i 
0.480058 - 0.075146i 
0.407175 - 0.064140i 
0.317805 - 0.050382i 
0.261841 - 0.041644i 
0.189994 - 0.030314i 
0.060091 - 0.009619i 


0.602 - 0.093Z 
0.545 - 0.085i 
0.481 - 0.076i 
0.408 - 0.064i 
0.318 - O.OSli 
0.2621 - 0.0417i 
0.1901 - 0.0303i 
0.06009 - 0.00962i 


0.599443 - 0.092703i 
0.543115 - 0.084496i 
0.480058 - 0.075146i 
0.407175 - 0.064140i 
0.317804 - 0.050382i 
0.261843 - 0.041643i 
0.189994 - 0.030314i 
0.060091 - 0.009619i 


0.602 - 0.093i 
0.545 - 0.085i 
0.481 - 0.076i 
0.408 - 0.064i 
0.318 - 0.051i 
0.2621 - 0.0417i 
0.1901 - 0.0303i 
0.06009 - 0.00962i 
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QNMs of gravitational field: I — 3, n — 1 



A 


u (WKB) (odd) 


CO (P-T) (odd) 


u (WKB) (even) 


u (P-T) (even) 




0.02 
0.04 
0.06 
0.08 
0.09 
0.10 
0.11 


0.58264 - 0.28129i 
0.53074 - 0.25536i 
0.47166 - 0.22639i 
0.40217- 0.19281i 
0.31550 - 0.15125i 
0.26056 - 0.12498* 
0.18952 - 0.09095i 
0.06008 - 0.02886i 


0.602 - 0.280i 
0.545 - 0.255i 
0.481 - 0.227i 
0.408 - 0.193i 
0.318 - 0.152i 
0.2621 - 0.1251i 
0.1901 - 0.0910i 
0.06009 - 0.02886i 


0.58264 - 0.28129i 
0.53074 - 0.25536i 
0.47166 - 0.22639i 
0.40217- 0.19281i 
0.31549 - 0.15125i 
0.26057 - 0.12497* 
0.18952 - 0.09095* 
0.06008 - 0.02886i 


0.602 - 0.280i 
0.545 - 0.255i 
0.481 - 0.227i 
0.408 - 0.193i 
0.318 - 0.152i 
0.2621 - 0.1251* 
0.1901 - 0.0910* 
0.06009 - 0.02886* 



QNMs of masslcss Dirac field: k = 1, n = 



A 


u (WKB) (-) 


co{ 


P-T) (-) 


u (WKB) (+) 


CO (P-T) (+) 





0.183 - 0.097* 


0.183 


- 0.098* 


0.183 - 0.095i 


0.189 - 0.105i 


0.02 


0.167- 0.087* 


0.166 


- 0.088* 


0.167- 0.085* 


0.171 - 0.094i 


0.04 


0.149 - 0.077* 


0.147 


- 0.077* 


0.149 - 0.075* 


0.151 - 0.082* 


0.06 


0.128 - 0.069i 


0.127 


- 0.065* 


0.128 - 0.064i 


0.128 - 0.069* 


0.08 


0.1005- 0.0507* 


0.0985 - 


- 0.0505i 


0.1006 - 0.0501i 


0.1001 - 0.0529i 


0.09 


0.0832 - 0.0418i 


0.0813- 


-0.0415i 


0.0833 - 0.0414i 


0.0824 - 0.0432i 


0.10 


0.0606 - 0.0304i 


0.0592 - 


-0.0301* 


0.0606 - 0.0303* 


0.0597- 0.0311* 


0.11 


0.0192 - 0.0096* 


0.0188- 


- 0.0096* 


0.0192 - 0.0096* 


0.0189 - 0.0097* 



QNMs of massless Dirac field: k = 2, n = 



A 


u (WKB) (-) 


u (P-T) (-) 


u (WKB) (+) 


u (P-T) (+) 




0.02 
0.04 
0.06 
0.08 
0.09 
0.10 
0.11 


0.38004 - 0.09642i 
0.34491 - 0.08715i 
0.30542 - 0.07690* 
0.25953 - 0.06516* 
0.20296 - 0.05085* 
0.16738 - 0.04190i 
0.12157- 0.03041i 
0.03849 - 0.00962i 


0.382 - 0.097i 
0.346 - 0.088* 
0.306 - 0.077* 
0.260 - 0.065* 
0.203 - 0.051* 
0.167 - 0.0418i 
0.121 - 0.0304* 
0.0384 - 0.0096* 


0.38007- 0.09637i 
0.34493 - 0.08712i 
0.30543 - 0.07688* 
0.25954 - 0.06515* 
0.20296 - 0.05084* 
0.16738 - 0.04190i 
0.12157- 0.03041i 
0.03849 - 0.00962i 


0.386- 0.099i 
0.349 - 0.089i 
0.308- 0.079* 
0.261 - 0.066* 
0.204 - 0.052* 
0.168 - 0.0423* 
0.122 - 0.0305i 
0.0384 - 0.0096* 



The relative error of the WKB method depends significantly on a field under consideration. 
Comparison of the sixth and fifth WKB order let us to conclude that the worst WKB convergence 
gives the Dirac field, where the expected error reaches 1% in the sixth WKB order. The results 
calculated in [10] by using the third WKB order differ from the sixth order values up to 4%. Thus, 
the possible relative error of the third order WKB results may exceed 4%. 

We can also evaluate the relative error order for gravitational modes by comparison QNMs of 
odd and even perturbations. Because the potentials can be expressed in the form 

VoM = W^ + ^+(3; V;.en = W;-^ + P; (13) 
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where 



2M 3 + 2c 3c2 + 2c3-9AM2 1 SAM^, ^ c^fc+l)^ 
= ^ + o , ^ T777(c +c ); 



r2 3r 3c(3M + cr) 3M' c " QM^ ' 

any solution for even parity can be found from the solution for odd parity [17]: 

*e^en(r*) = q (^W^g(r*) - ^odd(r*); q = const. 

Obviously, if a solution '^odd describing quasi-normal oscillations satisfies the boundary conditions 
(4), the '^even solutiou satisfics the same ones describing quasi-normal oscillations too. 
The same fact [12] takes place for massless Dirac QNMs of opposite chiralities: 



V.-K±^; WAr)^'-^. (14) 

Therefore, the QNM spectrum is the same for odd and even gravitational perturbations and for 
opposite chirality of Dirac perturbations in SdS background as well as in asymptotically fiat universe. 
We sec that the difference between the numerical results does not exceed our expectations. 

We can sec that the imaginary parts of QNMs, calculated using both sixth-order WKB and 
Poshl-Teller approaches, depend on multipole and overtone number and black hole mass only near 
the extremal value of A. This agrees with [11], where the QNMs of near extremal SdS space-time 
were found: 

uh^- (n + I + ^/(Z + 1)-^ (15) 
for scalar and electromagnetic perturbations. 



a;6=-(n + ^)z + J(Z + 2)(Z-l)-^ (16) 



for gravitational perturbations. Here 



h = —r-, — r, Tr 3M. 



(re-3M)(re + 6M)' 
Using the same technique, one can find that 



uh^-{n, + ]^^+^jK^-\ (17) 

for massless Dirac field perturbations. It means in the extremal limit fields of different spin decay 
with the same rate. Note, that formulae (15), (16), (17) give right results for any overtone number 
n and may be useful for studying the n — > oo limit. 

The well-known approximate formula for large multipole number I (or k in (9)) [18, 19]: 

^ = ^7^Al+\-{^+\)^)+0{\)- (18) 



3V3M 
1 

3V3M V V ' 2; 7 ' ^ \k 



u = — ^ — (k- (n+^i \ +0 (- 
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can be generalized for the case of non-zero A. Making use of the first-order WKB method or the 
formula of the approximation by the Poshl-Teller potential one can find for the non-extremal SdS 
black hole: 



VI - 9M2A /, 1 / 1\ \ , ^ 

u = j= (l + i:- [n + -]i]+0[-]] (19) 



3V3M V 2 V 2) J \l 
VI - 9M2A / / 1\ .\ ^/l 



4 Conclusion 

The QNMs of Schwarzschild-de Sitter black holes for fields of different spin have been calculated. 

The modes are determined by the black hole mass M and the cosmological constant A only. The 
frequencies all have a negative imaginary part, which means that the black hole is stable against these 
perturbations. The presence of the cosmological constant leads to decrease of the real oscillation 
frequency and to a slower decay. The interesting problem that was outside our consideration is 
the search of asymptotic behavior of SdS QNMs at large imaginary part {n — > oo) for which other 
approaches should be explored [20] . 
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